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ABSTRACT 

The ~70 unidentified sourees of the EGRET sky survey may be one of its most important legaeies. 
The identifieation of these sourees at other wavelengths is eritieal to understanding their nature. Many 
have flat speetra out to 10 GeV whieh, if extrapolated to TeV energies, would be easily deteetable 
relative to the steeply falling diffuse baekground. The Whipple Observatory y-ray teleseope has been 
used to observe a number of these whieh were seleeted based on their position, intensity and speetrum 
and in some eases based on a possible assoeiation with a supernova remnant or pulsar. No signifieant 
emission has been deteeted from these sourees, and upper limits are given. 

INTRODUCTION 

Despite extensive searehes for eounterparts to the ~30 low-lattitude unidentified EGRET sourees, 
the nature of these objeets is still largely unknown. Kaaret & Cottam (1996) have suggested that the 
low latitude unidentified sourees show a eorrelation with OB assoeiations, the sites of massive star 
formation. Sinee nearly half of all supemovae oeeur as the eore eollapse of young massive stars whieh 
explode into star formation regions (e.g., Huang & Thaddeus 1986) a eorrelation with the positions 
of pulsars and with supernova remnants also follows. Kaaret & Cottam argue that pulsars emit y- 
rays over a signifieantly longer lifetime than SNR, so that the number of y-ray pulsars is expeeted 
to be signifieantly larger than the number of y-ray SNR. However, Esposito et al. (1996), Stumer & 
Dermer (1994) and Stumer, Dermer and Mattox (1996) have presented evidenee for assoeiations of a 
number of these objeets with SNRs (y-Cygni, IC443, W44, Monoeeros) for whieh there is no pulsar 
within the 95% eonfidenee error eontour (Stumer, Dermer & Mattox 1996). Attempts to deteet radio 
pulsars in the error boxes of EGRET unidentified sourees have been unsueeessful (e.g., Niee & Sayer 
1997) and provide some eonstraints on models in whieh all of the Galaetie unidentified sourees are 
pulsars. Sinee EGRET generally laeks the spatial resolution to distinguish the point-like emission 
from pulsars and AGNs from the extended emission expeeted to arise in the vieinity of supernova 
shells, variability has been used to distinguish eompaet sourees. Dramatie transient sourees sueh as 
the enigmatie 2CG 135-1-1 and newly diseovered GRO J1838-I-04 are diffieult to interpret as either 
arising from AGNs or from pulsars, and are possibly representatives of a new elass of Galaetie y-ray 
souree distinet from isolated pulsars (e.g., Tavani et al. 1997). 


identification of these objects in two important ways. First, the contribution from the diffuse y-ray 
background falls as a steeper power of energy (~ than the source spectrum for many 

of these objects, implying a smaller effect from uncertainties in the diffuse background model in 
determining the position, flux and spectra of these sources. Second, the 0.13° angular resolution of 
the Whipple 10m y-ray telescope (Lessard and Buckley 1997) provides the ability to resolve extended 
sources such as SNRs and offers the potential to narrow the error box for bright sources. 

OBSERVATIONS AND ANALYSIS 

The high-energy y-ray telescope (Cawley et al. 1990) at the Whipple Observatory employs a 10 m 
diameter optical reflector to image Cerenkov light from air showers onto an array of 109 fast pho¬ 
tomultipliers covering a 3° field of view (FOV). By making use of the distinctive differences in the 
angular distribution of light and orientation of the shower images a y-ray signal can be extracted from 
the large background of hadronic showers. 

Data are generally taken in a differential mode where each 28 min ON-source run is followed 
by a 28 min OFF-source control run which is offset in right ascension to ensure that the same range 
of azimuth and zenith angles are sampled. While this cancels the zenith angle dependence of the 
cosmic-ray rate as well as other systematic effects in the camera, differences in sky brightness between 
the ON-source and OFF-source regions can lead to biases. For some galactic plane sources, such 
differences are substantial due to either diffuse emission from the galactic plane or bright stars. 

Systematic effects arising from such brightness differences can be largely canceled by the pro¬ 
cedure of software padding (Cawley et al. 1983). This procedure consists of adding noise to all pixels 
of each event so that matching PMTs in ON-source and OFF-source runs have identical noise pulse 
height spectra. Only PMT signals which exceed some multiple of this noise level are included in the 
subsequent analysis of the shower images (Punch et al. 1993). 

The technique used to generate the two-dimensional plots and upper limits is a simple extension 
of the standard Whipple data analysis (Reynolds et al. 1993) and is described in more detail in Lessard 
and Buckley (1997). After initial processing of the shower images including pedestal subtraction, 
gain correction and image cleaning (e.g.. Punch et al. 1993) individual Cerenkov shower images are 
subjected to a moment analysis to determine a set of parameters that characterize the roughly elliptical 
images. Each point on a two dimensional grid covering the 3° FOV is considered as the potential 
source position. For each event, the RMS width and length, centroid position, orientation, ellipticity 
and the skew of the shower image are calculated about this point of origin and tested for consistency 
with the parameter values expected for a y-ray event coming from the corresponding direction in the 
sky. For each grid point the number of candidate events ON-source and OFF-source are calculated, 
and the significance of the excess, S^, is derived using the likelihood ratio method of Li & Ma (1983). 
In the two-dimensional plots presented in Figures 1, the gray-scale indicates the number of excess 
counts (candidate y-rays) consistent with each grid point and the contours shown correspond to the 
statistic Sx in steps of 1. Note that due to the large number of trials associated with the 30x30 grid, 
approximately one Sx='i excess is expected for each two-dimensional plot. 

While it is desirable to have one control (OFF-source) run for each run ON-source, for some 
of the data presented here the number of exposures taken OFF-source is less than the number of 
ON-source runs. In this case, the background level is determined by normalizing the OFF-source 
data to the ON-source data in a 0.25° band around the perimeter of the field of view. The resulting 
normalization factor a enters into the calculation of the significance and the upper limit following the 
procedure of Li and Ma (1983). For the sources J0542-I-26, J0635-I-0521, and J1825-1307, little or no 
OFF-source data was taken and a background template was formed using contemporaneous control 
data taken for other sources. This results in an additional systematic error for these sources. 

In calculating upper limits, we are testing the hypothesis that the emission is coming from a 
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point source within the EGRET error box. Upper limits ealeulated for the extended regions eorre- 
sponding to the IC443, W44 and y-Cygni SNRs (for 2EG J0618+2234, 2EG J1857+0118, and 2EG 
J2020+4026 respeetively) are reported elsewhere (Buekley et al. 1997). 99.9% eonfidenee-level upper 
limits are ealeulated for eaeh grid point lying within the EGRET 95% eonfidenee interval using the 
method of Helene (1983) and aeeounting for the deelining y-ray deteetion effieieney away from the 
eamera eenter (Eessard & Buekley 1997). The maximum upper limits for eaeh error box are shown 
in Table 1. 

RESUETS AND CONCEUSIONS 

Table 1 shows preliminary Whipple upper limits for a number of unidentified sourees together with 
the extrapolated EGRET flux derived from the EGRET speetrum. Eluxes and speetral indiees are from 
Eierro et al. (1996). In addition to sourees from the first (Eiehtel et al. 1994) or seeond EGRET eata- 
logs (Thompson et al. 1995), we also inelude the souree J0749+17 from the initial list of unidentified 
sourees by Hartman et al. (1992). This objeet was not ineluded in the first EGRET eatalog beeause of 
its low signifieanee (< 4a), but is of interest sinee it prompted the radio pulsar seareh by Eundgren, 
Zepka and Cordes (1995) that led to the diseovery of the binary milliseeond pulsar PSR 0751+1807. 
Also on our list is the souree J0542+26 whieh was on the list of high eonfidenee unidentified sourees 
in the first but not the seeond EGRET eatalog. This souree has a 158 aremin error radius whieh eould 
not be shown in Eigure 1. This souree is of interest sinee it is eoineident with the position of the old, 
nearby (0.8-1.4kpe, Kundu et al. 1980) SNR S 147 as pointed out by Stumer and Dermer (1994). 

These data were taken over the period Deeember 1994 to May 1997. Two-dimensional plots for 
these sourees are shown in Eigure 1 exeluding results for J0618+2234, J1857+0118 and J2020+4026 
whieh are shown elsewhere (Buekley et al. 1997). Upper limits are at energies above 400 GeV unless 
otherwise indieated. 2EG 1746-2852 transits at an elevation of <30° resulting in an inerease of the 
effeetive area and energy threshold by a faetor of approximately 5.0 eompared with observations at the 
zenith (Krennrieh et al. 1997). While 2EG J1746-2852 shows a small (2.5o) exeess at the position of 
Sgr A* and within the EGRET error box, this exeess is not eonsidered signifieant given the additional 
systematie errors present for galaetie plane sourees. 2EG J0241+6119 shows a similar exeess near the 
position of 2CG135+01 and within the EGRET error box. The exeess in J0542+26 lies outside and to 
the south of the radio shell of S147 and approximately 0.5° away from the X-ray binary 4U0535+262, 
too far to make an assoeiation. The other sourees show no signifieant emission within the EGRET 
error boxes. Eurther deep observations with the GRANITE-III high resolution eamera should provide 
better sensitivity given the extended EOV and finer pixelization, and eorrelations with data taken at 
other wavelengths should improve ehanees for deteeting variable unidentified sourees. 
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Fig. 1: Plots of two-dimensional distributions of candidate gamma-ray events and contours of Sx 
for (a) J0241-\-6119 (cross at the position of 2CG135-\-01), (b) J0542-v26, (c) J0545-\-3943, (d) 
J0635-V0521, (e) J1746-2852 (cross at Sgr A*) and (f) J1825-v2234 (cross at the position of PSR 
B1823-13). Dotted contours show the elliptical fit to the EGRET 95% confidence intervals. . 
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